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Steady state performances of the electric power distribution system are 
normally assessed or evaluated based on load flow analysis. To properly carry 
out the analysis, a valid steady state load flow model of each distribution 
system component, including the wind power plant (WPP), needs to be 
developed. The present paper proposes a method for modeling and integrating 
squirrel cage induction generator (SCIG)-based variable speed WPP into a 
three-phase unbalanced distribution system load flow (DSLF) analysis. The 
proposed method is based on a single-phase T-circuit model of fixed speed 
WPP, which has successfully been applied to balanced electric power systems. 
In the present work, the single-phase T-circuit model is extended and modified 
to be used in steady state load flow analysis of three-phase unbalanced 
distribution systems embedded with SCIG-based variable speed WPP. Results 
of the case studies confirm the validity of the proposed method. 
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1. INTRODUCTION 

Load or power flow analysis is probably the most fundamental analysis in electric power system 
studies. Based on the analysis, the steady state performance of the electric power system can be assessed or 
evaluated. It has been acknowledged that a valid steady state model of each power system component needs to 
be derived and developed to enable the analysis to correctly be carried out. In the context of wind power plant 
(WPP) steady state load flow modeling, many methods have been proposed and developed in the last few years. 
Some of the methods are focused on fixed or near constant speed WPP modeling, and some deal with variable 
speed WPP modeling. 

Several techniques to incorporate fixed speed WPP have been introduced in some publications. The 
WPP steady state models in [1]—[3] are derived based on the equivalent circuit of the WPP induction generator. 
Three-node model of an asynchronous generator for load flow analysis has been proposed. The three-node 
model transforms the asynchronous generator equivalent circuit to several buses and branches [1]-[3]. A 
standard load flow program is then employed to obtain the load flow solution. A PQ model for an induction 
generator to be used in steady state load flow analysis has been developed [4]-[6]. The mathematical model in 
[4]-[6] is derived based on a formula that calculates the active and reactive powers of the induction generator. 
In the calculation of the powers, the turbine's mechanical power is assumed to be known. Various fixed speed 
WPP models for steady state load flow analysis have been introduced in [7], [8]. The WPP is modeled under 
unbalanced system conditions [8]. The developed model in [7] has the power formulations in terms of induction 
machine parameters, system terminal voltages, and mechanical input power. A method to integrate WPP into 
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a three-phase load flow analysis has been proposed in [8]. The method in [8] is developed by modifying the 
single-phase two-port network model so that it can be applied in three-phase systems. 

The methods discussed in [1]-[8] can only be applied to fixed speed WPP. Methods for modeling and 
integrating variable speed WPP have also been proposed and reported in many publications. A load flow model 
for doubly fed induction generators (DFIGs) is developed [9]. The electric power system and DFIG constraints 
are simultaneously solved in the model development. Also, in the method [9], the DFIGs are treated as PQ buses, 
and load flow analysis is carried out to obtain the stator voltages. An iterative technique for integrating DFIG- 
based WPP for load flow analysis has been introduced in [10]. The technique is based on DFIG steady state 
equivalent circuit. Moreover, various limits related to the DFIG are also taken into account in the method 
discussed in [10]. A scheme for incorporating the DFIG-based WPP in distribution systems load flow analysis 
has been proposed [11]. The proposed model is based on the DFIG equivalent circuit. The DFIGs are treated as 
PQ buses, and the DFIG reactive power limits are also taken into account in the method [11]. The forward- 
backward sweep (FBS) technique is then used to obtain the load flow solution of the power system networks [11]. 

The DFIG model for unbalanced distribution system load flow analysis has been proposed [12]. The 
model in [12] is derived based on an analytical representation of the WPP major components. The proposed model 
in [13] is also developed for three-phase load flow analysis. Effects of different rotor speeds are also considered 
in [13]. Steady state performance analysis of power flow control in DFIG is presented in [14]. Three different 
rotor speed modes of operation have been included in the analysis [14]. Simple techniques for modeling DFIG- 
based WPP for power system load flow studies have been proposed in [15]-[17]. In the method [15]-[17], the 
WPP model is developed based on formulas that calculate the WPP powers. Modeling of a power electronic 
converter is not needed in the method [15]-[17]. The steady state load flow model of SCIG-based variable speed 
WPP has been discussed [18]. The model in [18] is also derived based on formulas that calculate the WPP powers. 
On the other hand, in [19], a model for steady state and dynamic study of a grid-connected SCIG-based WPP has 
been presented. The model in [19] is established and implemented using MATLAB software. 

In most of the methods proposed in [1]-[19], the electric power systems have been assumed to be 
balanced. As the systems are balanced, single-phase approaches can be used to analyze the systems. It has been 
acknowledged that power systems are inherently unbalanced. Neglecting the unbalanced, especially if it is of 
a high level, will lead to invalid and unacceptable results. For such a system, three-phase methods of analysis 
should be used instead to obtain accurate results. Steady state load flow WPP models for three-phase 
unbalanced power systems have been introduced [7], [8]. However, the models are only valid for fixed speed 
WPP and cannot be applied to variable speed WPP. Therefore, against the above background, the present paper 
proposes a method for modeling and integrating SCIG-based variable speed WPP into a three-phase 
unbalanced DSLF analysis. The proposed method is based on a single-phase T-circuit model of fixed speed 
WPP, which has successfully been applied to balanced electric power systems [18]. In the present work, the 
single-phase T-circuit model is extended and modified to be used in steady state load flow analysis of three- 
phase unbalanced distribution systems embedded with SCIG-based variable speed WPP. Some case studies are 
also presented in this paper. In the case study, the validity of the proposed method is investigated and 
confirmed. 


2. THREE-PHASE DSLF PROBLEM FORMULATION 
In terms of nodal variables, the formulation of three-phase DSLF problem is given by (1) [8]. 


{[diag(V%*)]-1(sa¢ —_ Savery 23 yobcy abe =0 (1) 


abe 


where Vis nodal voltage vector, Sc” is power injection (generation) vector, S;“”° is power demand vector 
(load), and Y“"* is power system nodal admittance matrix. In (1), V, So’, Si@°, and Y” are of the forms: 


yabe = [yabe yabe .., yabeyT (2) 

sgh = (sage sage Sahel” @) 
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where n is the number of power system nodes, and: 


yore = (ye Ve eae (6) 
SH° =(S& SG SEIT (7) 
Se =(Sh Sh Stl? (8) 


yge y,eb yee 
uL u 
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In conventional DSLF analysis, power system nodes can be classified into two categories, namely: 
substation and load nodes. In the analysis, the substation node is usually considered the reference node. 
Therefore, this node's voltage magnitude and angle are specified at some values. The substation voltage 
magnitude is generally set at 1.0 pu (or a value near 1.0 pu), and the substation voltage angle is assumed to be 
zero (reference node). Since the distribution system is typically fed from the substation, the remaining nodes 
in the system can be considered load nodes. Therefore, the power generation of these nodes is all zeroes (see 
Table 1). 


Table 1. Known and unknown quantities 


Node Known quantity Unknown quantity 
Substation 1.020° sge 
ve" = 11.02 —120° 
1.02120° 
Load Sane = i) 0 0]? pare 


3. MODEL OF SCIG-BASED VARIABLE SPEED WIND TURBINE 

It has been acknowledged that variable speed wind turbines can harness wind energy more efficiently 
than fixed speed wind turbines. This advantage is probably the main reason why variable speed wind turbines 
have been replacing fixed speed wind turbines in the last few decades. Several configurations can be used to 
implement the variable speed operation of a wind turbine. However, due to its relatively lower cost and easy 
maintenance, the SCIG-based variable speed wind turbine is gaining popularity nowadays. 

Figure 1 shows a basic structure of a SCIG-based variable speed wind turbine system [18]-[23]. It 
can be seen from Figure | that the main component of the system is SCIG which converts turbine mechanical 
power to electrical power. The system also contains a power electronic converter (PEC), which consists of a 
power converter/inverter, and a DC-link. The PEC will isolate the SCIG speed from the frequency of the power 
grid where the WPP is connected. This isolation makes the WPP operation speeds wider, and the wind energy 
can be extracted more efficiently. Moreover, with the presence of PEC, the WPP can be set to control the 
reactive power output and support the system's reactive power demand. 

The equivalent circuit of SCIG is shown in Figure 2 [18]-[23]. In Figure 2, the electrical quantities 
Vs, Is, Vr, and Ip denote stator voltage, stator current, rotor voltage, and rotor current, respectively. The circuit 
impedances Zs, Zr, and Zy in the figure are determined using: 


Z5 = Rs + jXs (11) 
Zp = Rat jXr (12) 
Zum =IRXm/(Re + iXm) (13) 


where Rs, Xs, Rr, Xr, R., and X,, are the stator circuit resistance, stator circuit reactance, rotor circuit resistance, 
rotor circuit reactance, magnetic core circuit resistance, and magnetic core circuit reactance, respectively. 
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Figure 2. Steady state equivalent circuit of SCIG 
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In [18], it has been shown that a single-phase steady state model of SCIG-based variable speed WPP 


valid for a load flow analysis of a balanced power system is of the form: 


Zri1Vy + Zr2Ss + Zr355 + Zr4Sy — Py = 0 


where: 
x 
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(14) 


(5) 


(16) 


(17) 


(18) 


(19) 


(20) 


In (15) and (16), 7. is the PEC efficiency. It is to be noted that the power injection (i.e., power 
generation) at the WPP bus is also the WPP power output. In the present work, it will be assumed that the WPP 
is operated at a unity power factor, i.e., there is no reactive power delivered/absorbed to/from the power system 


(Q, = 0). Therefore, power injection at the WPP bus is: 


Se = Fy = 1cPs 


(21) 


4. THREE-PHASE MODEL OF SCIG-BASED VARIABLE SPEED WIND TURBINE 
To facilitate the load flow analysis in a three-phase system, the WPP single-phase model (14) is 


extended to become: 
Zabeyare 4 zabecabe 4 zabe(gabcy: 4 Zabegabe _ pabe — 9 
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where: 
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In (27), the Z7“” sub-matrices are determined based on (17) - (20) as: 
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where p = a, b, c. In (28) - (31), Zs’, Zr’, and Zy’ are calculated based on (11) - (13) as: 
P 

Zs = RE + jXS (32) 

Zp = Re + jXp (33) 

Zim = IREXm/ (Re + iXm) (34) 


According to (1) and (22) represent the complete formulation of the three-phase load flow problem of 
the electric power distribution system embedded with SCIG-based variable speed WPP. The quantities to be 
calculated and the equations to be solved are given in Table 2. By looking at Figure 1, it is clear that at the 
WPP node, the following relationships are valid: 


yore = ye (35) 
SGPo = Sq = NPP" (36) 
where: 
T 
Vee — [V," VP Vie (37) 
aa ae (38) 
pyc = [P$ PS PS]? (39) 
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Table 2. Known and unknown quantities for system with WPP 


Node Equation Known quantity Unknown quantity 
Substation (1) 1.020° ser" 
v7" = |1.02 — 120° 
1.02120° 
Load (1) sz =[0 0 oj? yare 
WPP (1) & (22) Pre’; 2832 Ri Ln ye" soe 


5. CASE STUDY 
5.1. Case study 1 
5.1.1. Test system 

To confirm the method proposed in section 4, a 12-node balanced distribution system 
shown in Figure 3 will be used in this first case study [24], [25]. Although the system is balanced and can be 
analysed by single-phase approach, for the purpose of validation study, it will be analysed using the proposed 
three-phase method. The system voltage in Figure 3 is 11 kV with a total three-phase load of 4.305 kW and 
1.215 kKVAR. The system data is presented in Table 3. In this case study, WPP with a power rating of 3.5 MW 
is assumed to be connected to node 12 via a step-up transformer. Data of the WPP is given in Table 4. All data 
in pu are on 1 MVA base. 


5.1.2. Turbine mechanical power 

It has been acknowledged that turbine mechanical power (P,,,) depends highly on the wind speed. 
Table 5 shows the amount of turbine mechanical powers for various wind speeds used in the first case study. 
The turbine powers have been calculated using the formula given in [15]. In the calculation, the values of air 
density and turbine performance coefficient are assumed to be 1.225 kg/m? and 0.45, respectively. 


5.1.3. Results and discussion 

Tables 6-9 present the results of the load flow analysis for the 12-node distribution system. Because 
the system is balanced, electric quantities in the three phases are all identical. It can be seen that the results in 
Tables 6-9 are in excellent agreement with the results in [18] which confirm the proposed method validity. 


Substation 


WPP 


Figure 3. Test system (12-node network) 


Table 3. Test system data (12-node network) 
Line No. Sending Bus _ Receiving Bus R* x* Per Oe 
1 1 2 1.093 0.455 60 60 


2 2 3 1.184 0.494 140 30 
3 3 4 2.095 0.873 155 55 
4 4 5 3.188 1.329 130 30 
5 5 6 1.093 0455 120 15 
6 6 7 1.002 0417 155 55 
7 7 8 4.403 1.215 145 45 
8 8 9 5.642 1.597 140 40 
9 9 10 2.890 0.818 135 30 
10 10 i 1.514 0.428 140 30 
ul id 12 1.238 0.351 115 15 


*Line resistance/reactance in ohm/phase 
**Toad connected to receiving bus in kW and kVAR (perphase) 
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Table 4. WPP data (12-node network) 
Component Quantity 
Turbine Blade length: 45 meter 
Rated power: 3.5 MW 
Speed: 
Cut-in: 5 m/s; Rated: 12 m/s; Cut-out: 23 m/s 
Gear Box Ratio: 1/60 
Generator Type: SCIG 
Rated power: 3.5 MW 
Pole pairs: 4 
Voltage: 690 Volt 
Resistances/Reactances (in pu): 
Ro = RP = Rf = R= RP = Rf = 0.01 
Xf=X) =X6=X2=X! =X¢ = 0.10 
R= RE = RE = 500, Xp! = Xu? = Xp = 50 
PEC Efficiency: 95% 
Step-Up Impedance (in pu*): Z’ = Z° = Z = j0.1 
Transformer 
Table 5. Turbine power (12-node network) 
Wind Speed (m/s) _ Turbine power (MW) 
5 0.3367 
6 0.5346 
7 0.7980 
8 1.3362 
9 1.5586 
10 2.0745 
11 2.6933 
12 3.4243 
Table 6. Power output of WPP (12-node network) 
Pn(MW) — Phase-a(MW) _ Phase-b (MW) _ Phase-c (MW) —_ 3-Phase (MW) 
0.3367 0.1035 0.1035 0.1035 0.3105 
0.5346 0.1640 0.1640 0.1640 0.4920 
0.7980 0.2435 0.2435 0.2435 0.7305 
1.3362 0.4023 0.4023 0.4023 1.2069 
1.5586 0.4668 0.4668 0.4668 1.4004 
2.0745 0.6140 0.6140 0.6140 1.8420 
2.6933 0.7868 0.7868 0.7868 2.3604 
3.4243 0.9862 0.9862 0.9862 2.9586 
Table 7. WPP terminal voltage magnitude (12-node network) 
P.»(MW)__Phase-a(pu) — Phase-b (pu) _ Phase-c (pu) 
0.3367 0.8902 0.8902 0.8902 
0.5346 0.9059 0.9059 0.9059 
0.7980 0.9254 0.9254 0.9254 
1.3362 0.9614 0.9614 0.9614 
1.5586 0.9751 0.9751 0.9751 
2.0745 1.0045 1.0045 1.0045 
2.6933 1.0363 1.0363 1.0363 
3.4243 1.0701 1.0701 1.0701 
Table 8. Distribution line loss (12-node network) 

P,,(MW)  Phase-a (MW, MVAR) _ Phase-b (MW, MVAR) __ Phase-c (MW, MVAR) _3-Phase (MW, MVAR) 
0.3367 0.1418+j0.0563 0.1418+j0.0563 0.1418+j0.0563 0.4254+j0.1689 
0.5346 0.1206+j0.0504 0.1206+j0.0504 0.1206+j0.0504 0.3618+j0.1512 
0.7980 0.0976+j0.0455 0.0976+j0.0455 0.0976+j0.0455 0.2928+j0.1365 
1.3362 0.0650+j0.0436 0.0650+j0.0436 0.0650+j0.0436 0.1950+j0.1308 
1.5586 0.0561+j0.0454 0.0561+j0.0454 0.0561+j0.0454 0.1683+j0.1362 
2.0745 0.0438+j0.0545 0.0438+j0.0545 0.0438+j0.0545 0.1314+j0.1635 
2.6933 0.0416+j0.0726 0.0416+j0.0726 0.0416+j0.0726 0.1248+j0.2178 
3.4243 0.0726+j0.1021 0.0726+j0.1021 0.0726+j0.1021 0.1578+j0.3063 
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Table 9. Substation power (12-node network) 


P,,(MW) _ Phase-a(MW,kVAR) _ Phase-b (MW, MVAR) __Phase-c (MW, MVAR) _3-Phase (MW, MVAR) 
0.3367 1.4733+)0.4613 1.4733+)0.4613 1.4733+j0.4613 4.4199-+j1.3839 
0.5346 1.3916+j0.4554 1.3916+j0.4554 1.3916+j0.4554 4.1748-+j1.3662 
0.7980 1.2891+)0.4505 1.2891+)0.4505 1.2891+4j0.4505 3.8673+1.3515 
1.3362 1.0977+j0.4486 1.0977+j0.4486 1.0977+j0.4486 3.29314 1.3458 
1.5586 1.0243+j0.4504 1.0243+j0.4504 1.0243+j0.4504 3.0729+1.3512 
2.0745 0.8648+)0.4595 0.8648+)0.4595 0.8648+)0.4595 2.5944+j1.3785 
2.6933 0.6897-+j0.4776 0.6897-+j0.4776 0.6897+)0.4776 2.069141.4328 
3.4243 0.5015+j0.5071 0.5015+j0.5071 0.5015+j0.5071 1.5045+j1.5213 


5.2. Case study 2 
5.2.1. Test system 

To further validate the method proposed in section 4, a 19-node unbalanced distribution system shown 
in Figure 4 will be used in this case study [24], [25]. The system voltage in Figure 4 is 11 kV with a total three- 
phase load of 515.94 kW and 177.27 kVAR (Phase-a: 176.33 kW and 61.23 kVAR; Phase-b: 166.34 kW and 
56.34 kKVAR; Phase-c: 173.27 kW and 59.70 kKVAR). The system data is presented in Tables 10 and 11. In this 
case study, WPP with a power rating of 300 kW is assumed to be connected to node 19 via a step-up 
transformer. Data of the WPP is given in Table 12. All data in pu are on | MVA base. 


5.2.2. Turbine mechanical power 

Table 13 shows the turbine mechanical powers for various wind speeds used in the second case study. 
The turbine powers have also been calculated using the formula given in [15]. In the calculation, the values of 
air density and turbine performance coefficient are assumed to be 1.225 kg/m? and 0.35, respectively. 


Figure 4. Test system (19-node network) 


Table 10. Test system line data (19-node network) Table 11. Test system load data (19-node network) 


No. Line Self-impedance Mutual impedance Node Phase-a Phase-b Phase-c 
(Ohm) (Ohm) P (kW) Q(kVAR) P (kW) Q (kVAR) P (kW) QO (kVAR) 
1 1-2  3.0x(1.56090+4j0.67155) 3.0x(0.52030+4j0.22385) 1 0 0 0 0 0 0 
2 2-3  5.0x(1.560904j0.67155) 5.0x(0.52030+4j0.22385) 2 1038 5.01 5.19 2.52 10.38 = 5.01 
3 2-4 — 1.5x(1.56090+4j0.67155) 1.5x(0.520304j0.22385) 3 11.01 5.34 5.19 2.52 9.72 4.71 
4 4-5  1.5x(1.56090+4j0.67155) 1.5x(0.52030+j0.22385) 4 405 1.95 5.67 2.76 648 3.15 
5 4-6  1.0x(1.560904j0.67155) 1.0x(0.520304j0.22385) 5 648 3.15 5.19 252 453 2.19 
6 6-7 2.0x(1.56090+j0.67155) 2.0x(0.520304j0.22385) © 420 204 3.09 150 291 141 
7 6-8 2.5x(1.56090+j0.67155) 2.5x(0.52030+j0.22385) : ole es : - ue na se 
8 8-9 —3.0x(1.56090+4j0.67155) 3.0x(0.52030+j0.22385) 9 4539 507 1491 723 1329 642 
9 9-10 5.0x(1.56090+j0.67155) 5.0x(0.52030+)0.22385) 49 3,39 165 420 204 2.58 126 
11 10-12 1.5x(1.56090+j0.67155) 1.5x(0.52030+j0.22385) 12 9.72 4.71 8.10 3.93 8.10 3.93 
12 11-13 5.0x(1.560904j0.67155) 5.0x(0.52030+j0.22385) 13. 438 2.13 5.34 2.58 6.48 3.15 
13. 11-14 1.0x(1.56090+j0.67155) 1.0x(0.52030+j0.22385) 14. 3.09 1.50 3.09 1.50 4.05 1.95 
14. 12-15 5.0x(1.560904j0.67155) 5.0x(0.52030+j0.22385) 15 1438 2.13 1486 2.34 16.96 3.36 
15 12-16 6.0x(1.56090+4j0.67155) 6.0x(0.520304j0.22385) 16 17.77 3.78 2038 5.01 17.77 3.78 
16 14-17 3.5x(1.560904j0.67155) 3.5x(0.52030+j0.22385) 17 1648 3.15 1486 2.34 1486 2.34 
17 14-18 4.0x(1.560904j0.67155) 4.0x(0.52030+4j0.22385) 18 15.34 2.58 15.34 2.58 15.52 2.67 
18 15-19 4.0x(1.560904j0.67155) 4.0x(0.52030+j0.22385) _19 18.76 4.23 20.05 4.86 17.14 3.45 
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Table 12. WPP data (19-node network) 


Component Quantity 
Turbine Blade length: 16 meter 
Rated power: 300 kW 
Speed: 
Cut-in: 5 m/s; Rated: 12 m/s; Cut-out: 23 m/s 
Gear Box Ratio: 1/40 
Generator Type: SCIG 


Rated power: 300 kW 
Pole pairs: 2 
Voltage: 690 Volt 
Resistances/Reactances (in pu): 
Ro= Ro =Ro = RA = R? = RK = 0.0020 
Xf=X! =X = X= X) =X, = 0.0100 
Ro= RL = RE = 500, Xn = Xn? = Xn = 25 
PEC Efficiency: 95% 
Step-Up Transformer Impedance (in pu*): Z’ = Z’ = Z° = j0.1 


Table 13. Turbine power (19-node network) 
Wind Speed (m/s) _ Turbine Power (kW) | Wind Speed (m/s) _ Turbine Power (kW) 


P) 21.5513 9 125.6873 
6 37.2407 10 172.4106 
7 59.1368 11 229.4785 
8 88.2742 12 297.9255 


5.2.3. Results and discussion 

Tables 14-17 present the results of the load flow analysis for the 19-node distribution system. Because 
the system is unbalanced, it can be seen that the electric quantities in the three phases are not the same. The 
results of the load flow analysis also show that with the increase in turbine mechanical power input, the WPP 
active power output also increases (see Table 14 and Figure 5). However, its value is slightly smaller than the 
mechanical power due to power losses in the WPP (1.e., losses in SCIG and PEC). 

With the installation of WPP, one of the advantages is that the system voltage profile can be improved 
(see Table 15 and Figure 6). This improvement is caused by the WPP power injection. It can be seen that the 
voltage profile gets better with the increase in WPP power injection (i.e., WPP output). The system voltage 
profile improvement will, in turn, reduce the distribution line power losses as the power injection increases 
(see Table 16 and Figure 7). Another important advantage obtained with the presence of WPP is that the power 
supplied by the substation can be reduced (see Table 17 and Figure 8). The reduction in substation power is 
possible because some of the system loads can be fed by the WPP. It is also to be observed that for every case 
of turbine power value, the substation power plus WPP power always equals total system load plus distribution 
line losses. It is to be noted that the line losses have been calculated based on the line impedances and currents. 
These results further verify the validity of the method proposed. 


Table 14. Power output of WPP (19-node network) Table 15. WPP terminal voltage magnitude 
Pe, Phase-a_ Phase-b_  Phase-c —_3-Phase (19-node network) 

(kW) (kW) (kW) (kW) (kW) P,,(kW)  Phase-a(pu) _ Phase-b (pu) _ Phase-c (pu) 
21.5513 5.01 5.02 5.02 15.05 21.5513 0.92691 0.92492 0.92565 
37.2407 9.97 9.97 9.97 29.91 37.2407 0.93075 0.92879 0.92951 
59.1368 16.88 16.89 16.88 50.65 59.1368 0.93604 0.93412 0.93482 
88.2742 26.08 26.08 26.08 78.24 88.2742 0.94298 0.94111 0.94179 
125.6873 37.89 37.89 37.89 113.67 125.6873 0.95171 0.94990 0.95056 
172.4106 52.63 52.64 52.64 157.91 172.4106 0.96237 0.96062 0.96125 
229.4785 70.65 70.65 70.65 211.95 229.4785 0.97502 0.97335 0.97394 
297.9255 92.24 92.25 92.25 276.74 297.9255 0.98971 0.98812 0.98868 


Table 16. Distribution line loss (19-node network) 
P,,(KW) _ Phase-a (kW, KVAR) _Phase-b (kW, KVAR) __ Phase-c (kW, KVAR) _3-Phase (kW, KVAR) 


21.5513 9.574)4.16 9.68-4)4.12 9.81-4)4.23 29.06+j12.51 

37.2407 8.85+)3.86 8.95+j3.82 9.08-+j3.92 26.88+j11.60 
59.1368 7.944j3.48 8.00+j3.44 8.14+4j3.54 24.08+j10.46 
88.2742 6.86-+)3.06 6.89-+)3.01 7.04+4j3.10 20.78+j9.17 
125.6873 5.70+j2.64 5.69-+j2.58 5.85+j2.67 17.244j7.89 
172.4106 4.584)2.31 4.534j2.23 4.714)2.31 13.834j6.85 
229.4785 3.70-4j2.15 3.60-+4j2.07 3.79-4j2.13 11.094j6.35 
297.9255 3.284j2.31 3.12-4]2.22 3.34-4j2.27 9.744j6.80 
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Table 17. Substation power (19-node network) 

Py Phase-a Phase-b Phase-c 3-Phase 
(kW) (kW,KVAR) (KW, kKVAR) (KW, KVAR) (kW, KVAR) 
21.5513 180.89+j65.39 171.00+j60.46 = 178.06+j63.93 529.94+j189.78 
37.2407 —-175.224j65.09 165.314j60.16 172.384j63.62 512.904j188.87 
59.1368  167.394j64.71 —_157.454j59.78  164.534j63.24 489.374j187.73 
88.2742 157.114+j64.29 — 147.144+j59.35  — 154.234+j62.80 458.48+j186.44 
125.6873 144.144+j63.87 134.134j58.92  141.234+j62.37  419.51+j185.16 
172.4106 =128.28+j63.54 = 118.234j58.57 = 125.34+j62.01 371.85+j184.12 
229.4785 109.384j63.38 99.284j58.41 — 106.414j61.83 315.084j183.62 
297.9255 87.364j63.54 _77.214)58.56 _84.364j61.97 __248.934j184.07 
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6. CONCLUSION 

A method for modeling and integrating SCIG-based variable speed WPP into three-phase unbalanced 
DSLF analysis has been proposed in this paper. The method is based on a single-phase T-circuit model of fixed 
speed WPP, which has successfully been applied to balanced electric power systems. In the present work, the 
single-phase T-circuit model is extended and modified so that can be applied in steady state load flow analysis of 
three-phase unbalanced distribution systems embedded with SCIG-based variable speed WPP. Case studies have 
also been presented in this paper. In the case studies, validity of the proposed method is investigated and 
confirmed. Results of the case studies also demonstrate the capability of the proposed method in representing the 
SCIG-based variable speed WPP for steady state load flow studies of the electric power distribution systems. 
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